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Abstract 

This note is a correction to a paper of Cortez, Peskin, Stockie & Varela [SIAM J. 
Appl. Math., 65(2):494-520, 2004], who studied the stabihty of a parametrically-forced, 
circular, elastic fiber immersed in an incompressible fluid in 2D, and showed the ex- 
istence of parametric resonance. The results were represented as plots that separate 
parameter space into regions where the solution is either stable or unstable. We uncov- 
ered two errors in the paper: the first was in the derivation of the eigenvalue problem, 
and the second was in the code to used to calculate the stability contours. 

An error in the derivation was found in Appendix B of yj, where the gradient of the 
Dirac delta function was improperly transformed from Cartesian to polar coordinates. This 
led to an error in the perturbation expansion of the immersed boundary forcing term. The 
correct transformation is |2] 

V6{x - = ^V[5(r - 1)6(6 - s)], 
which changes the second equation in Claim 1 to 

\ r J r 
where subscripts denote partial derivatives. This in turn modifies equations (6.1)-(6.4) in [T] 
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for the eigenvalue problem. For the case when (7 + in) 7^ 0, the equations are 
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while when (7 + in) = we have 



= ip\X'^ - itXI_^ + zrX^+i) + v{Xl - ^rX;„i + irXl^^), 
= v\K - + ^rX^+,) + ip\X: - ^rX^.^ + zrX;+,). 



Note that in each pair of equations, the correction only affects the second equation while 
the first remains unchanged. These results were derived with the Navier-Stokes equations 
written in stream function-vorticity variables. The derivation was repeated in terms of 
primitive variables, and identical results were obtained. 

The second error was in the Matlab code used to generate the stability plots, which 
affected the stability plots in three ways: 

• One solution mode was missing for the lowest wavenumber "subharmonic" case. 

• Spurious points were obtained in the region where the subharmonic mode was supposed 
to have been located (i.e., for small values of wavenumber p), as well as in the upper 
right reaches of the plot (for large values of both wavenumber and amplitude, see 
especially the rightmost plot in Figure 8.2). There was also a spurious line of points 
running just above the line r = 0.5. 

• The borders of the stability regions were in some cases noticeably deformed. 

All of these discrepancies can be seen by comparing the original Figure |A] (which repeats the 
stability plots for cases I and II from Figure 9.1 in [1]) with the new Figure |B] (which shows 
the corrected plots generated using the new equations). Both of these plots correspond to 
the following parameter values: 

Case I Case II 

K = 0.5 K = 0.04 

u = 0.004 z/ = 0.00056 

= 3.2 X 10"^ = 7.84 X 10-6 

Each point depicted in these plots represents one solution to the eigenvalue problem 
formulated in Eq. (8.1) of [1], and the stabihty boundaries are traced out by varying the 
angular wavenumber [p) of the circular fiber. The blue points denote harmonic modes and 
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K = 0.5, v = 0.004, (|) = 3.2e-05 



K = 0.04, v = 0.00056, i|) = 7.g4e-06 
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Fig. 9.1. Stability diagrams for Case I (k = 0.5, v = O.OQi) and Case II (k = 0.04, u = 0.00056;. 



Figure A: Original stability contours for cases I and II, reproduced from [T]. 
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Figure B: Corrected stability contours for cases I and II (compare to Figure [a| . 



the red points denote subharmonic modes. The only "physical" instabilities are those corre- 
sponding to integer values of wavenumber p that lie inside the stability "tongues" and satisfy 
^ r ^ ^. With this in mind, the lowest wavenumber unstable mode in case I corresponds 
to p = 2, while that in case II has p = 4. 

Other than the discrepancies noted above, the original stability plots from pQ are still 
qualitatively correct and the parametric resonances identified in the original paper are also 
true instabilities (as verified in numerical simulations). Nevertheless, we have performed 
new calculations with the corrected Matlab code that demonstrate the existence of resonant 
subharmonic modes that were not identified in the original paper. New resonant modes 
occur only when the missing subharmonic region overlaps with an integer value of p, which 
is not the case in the two examples plotted in Figures |A] and |B} The original stability plots 
from Figures 8.2 and 8.3 are reproduced in Figure O while our corrected plots are given 



in Figure D In all cases, the spurious points are eliminated and a new subharmonic mode 
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Fig. 8.2. A series of three plots showing the impact of changes in the stiffness on the stability 
boundaries in the p,T-plane, with viscosity v = 0.0002 and three different values of stiffness: k, = 0.02 
(left), 0.04 (middle), and 0.08 (right). 
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Fig. 8.3. A series of three plots showing the impact of changes in viscosity in the stability 
boundaries in the p, T-plane, with stiffness k = 0.04 and three different values of viscosity: v = 
0.00005 (left), 0.001 (middle), and 0.005 (right). 



Figure C: Original stability plots, reproduced from 



appears. 

To illustrate the presence of the missing unstable modes, we choose the following two 
sets of parameters listed below 



Case III 



Case IV 



0.02 
0.0002 

2 X 10"*^ 



K 
V 



0.08 
0.0002 

5 X 10-^ 



which correspond (respectively) to the left- and right-most plots in Figures 8.2-8.3. Numerical 
simulations of the full immersed boundary equations have been performed for Case III (with 
p = 3) and Case IV (with p = 2) to verify that these parametric resonances are actual 
instabilities. Following the approach used in [1], we initialize the fiber using the perturbed 
circular shape r = i?(l -|- 0.05 cos(p^)), where p is chosen equal to the resonant wavenumber. 
Figure [E] shows the amplitude of various p-modes in the two simulations, from which it is 
clear that the resonant mode is excited and its amplitude grows in time, as would be expected 
for the resonant case. The amplitude of the other modes remains small. 

Other than the presence of these missing lowest-wavenumber subharmonic modes, the 
effect on the stability contours is relatively small for the parameter values we considered. 



4 



K = 0.02, V = 0.0002, <|) = 2e-06 



K = 0.04, V = 0.0002, i|) = le-06 



K = 0.08, V = 0.0002, $ = 5e-07 




1 2 3 4 5 6 7 
wavenumber, p 




K = 0.04, V = 5e-05, 41 = 6.256-08 



1 2 3 4 5 6 7 
wavenumber, p 



K = 0.04, V = 0.001 , 4. = 2.5e-05 



J 



1 2 3 4 5 6 7 
wavenumber, p 



K = 0.04, V = 0.005, = 0.000625 



1 2 3 4 5 6 7 
wavenumber, p 



Figure D: Corrected stability plots (compare to Figure O. 
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;ure E: Amplitudes of the various p-modes from numerical simulations of Cases III 
= 3 unstable) and IV (right, p = 2 unstable). 
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Figure F: Comparison of the eigenvalue plots for the original (blue) and corrected (red) 
equations in Cases I (left) and II (right). 



Figure [F] draws a direct comparison between the stability regions for cases I and II obtained 
with the original equations and with the corrected equations. Clearly, the correction results 
in only a small deformation in the contours and thus has only minimal impact on the solution. 
All other cases considered in [1] yield similar results. 
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